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Q\ ■ Abstract 

p i' We perform a systematic study of radiative corrections to the masses of the Higgs bosons in 

the minimal supersymmetric standard model (MSSM) augmented by a single gauge singlet, 
' the so-caUed next-to-minimal supersymmetric standard model (NMSSM). Our method is 
based on the one-loop effective potential and includes effects of top quark, squark, Higgs 



X 

■ and Higgsino loops. We discuss the renormalisation group flows of Yukawa couplings 
and the upper bound on the lightest CP-even neutral Higgs boson mass as a function of 
the heavier stop mass and top mass. We then give a general discussion of Higgs boson 
phenomenology including radiative corrections. We survey as much of the parameter space 
of the Higgs sector of the NMSSM as is practicable, and analyse the full spectrum of Higgs 
masses and couplings in these regions of parameter space. Characteristic signatures of the 
NMSSM such as light charged bosons and weakly coupled neutral scalars are discussed, as 
are the relative sizes of the various radiative corrections. The MSSM is also discussed as 
a limiting case of the NMSSM for comparison. 
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1 Introduction 

The most widely studied extension to the standard model is supersymmetry (SUSY) [1], in 
which the particles of the theory are supplemented by the inclusion of their superpartners. 
SUSY requires the introduction of a non- minimal Higgs sector with at least two doublets [2] . 
An important question is how heavy the lightest neutral CP-even supersymmetric Higgs 
boson, h, can be within the framework of supersymmetric grand unified theories (SUSY 
GUTs). In SUSY GUTs all the Yukawa couplings are constrained to remain perturbative 
in the region Msusy ~ 1 TeV to Mqut ~ 10^^ GeV. This constraint provides a maximum 
value at low energies for those Yukawa couplings which are not asymptotically-free, and 
is obtained from the renormalisation group (RG) equations together with the boundary 
conditions that the couplings become non-perturbative at Mqut - the so-called "triviality 
limit". In the minimal supersymmetric standard model (MSSM) [1, 2], the triviality limits 
provide a useful bound on the top quark mass m^. The upper bound on the h mass, 
m/i, in the MSSM, including radiative corrections, has recently been the subject of much 
discussion [3, 4, 5, 6, 7]. 

However the MSSM is not the most general low energy manifestation of SUSY GUTs. 
It is possible that SUSY GUTs give rise to a low energy theory which contains an addi- 
tional gauge singlet field, the so called next-to-minimal supersymmetric standard model 
(NMSSM) [8, 9, 10]. In all supersymmetric models, including those with non-minimal 
Higgs sectors, there is an upper bound on the lightest CP-even scalar Higgs mass, which 
is generally larger for non-minimal models [11, 12]. Although the NMSSM involves a sin- 
gle gauge singlet, the estimates of the bound are applicable to a model with an arbitrary 
number of Higgs singlets, since the singlet fields may always be redefined so that only one 
of them couples to the doublets and only this field contributes to the upper two-by-two 
mass matrix which gives the bound. This argument cannot be applied to models which 
contain extra non-singlet degrees of freedom (although these usually give a lower value 
for the bound) or to the description of the whole spectrum in any model other than the 
NMSSM. 

There has recently been much interest in radiative corrections to Higgs boson masses 
in the NMSSM [13, 14, 15, 16, 17, 18, 19, 20]. The two approaches that have been 
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considered are the renormalisation group (RG) approach [16, 17, 18], and the one-loop 
effective potential approach [13, 14, 15, 19, 20]. In the RG approach one derives an effective 
low energy Higgs potential at a low energy scale ^ from SUSY matching conditions at a 
scale MsuSY- The squarks and other sparticles are usually assumed to be degenerate 
at MsuSYi and the logarithmic radiative corrections are efficiently summed by the RG 
equations for the case of either one [16, 17] or two [18] light Higgs doublets. This approach, 
although convenient, fails to pick up non-logarithmic corrections, and in general becomes 
complicated when a general squark spectrum is considered [7]. The one-loop effective 
potential [21, 22] has been used in the MSSM to estimate the radiative corrections due to 
a general squark spectrum [4] and these calculations have been repeated in the NMSSM 
[14, 15, 19]. However in the NMSSM there are other corrections due to Higgs loops which 
may become large due to the presence of large couplings and trilinear soft parameters, and 
it is one of the purposes of the present paper to consider such contributions to the one-loop 
effective potential. 

In this paper we shall perform a systematic study of radiative corrections to Higgs 
boson masses in the NMSSM, using the one-loop effective potential. We shall consider 
the effects of all particles which couple through relatively large (order one) couplings in 
the model including the effects of loops of top quarks, stop squarks, Higgs bosons, and 
Higgsinos. The top and stop corrections have been calculated before [14, 15, 19] and we 
include them in our analysis for completeness. The Higgs and Higgsino contributions to the 
effective potential have not been calculated before, and are dealt with here by numerical 
techniques since their contribution has no simple analytic form. 

Having developed the techniques for dealing with radiative corrections in the NMSSM 
we then apply these techniques in two different ways. The first application is to the problem 
of the bound on the lightest CP-even Higgs boson mass, which was discussed above. More 
generally we give a phenomenological discussion of Higgs boson masses and couplings over 
as much of the parameter space of the NMSSM as it is feasible to consider. Our goal is to 
understand the behaviour of the Higgs boson spectrum as each of the parameters in turn is 
varied, and to give some examples of phenomenological signatures which would enable the 
NMSSM to be distinguished from the MSSM. It is important to stress that we shall not 
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impose GUT scale constraints on the soft SUSY-breaking parameters, so that our analysis 
is of a general low energy phenomenological nature. 

We shall also present an RG analysis of the dimensionless Yukawa couplings of the 
model which will turn out to provide a useful guide to the typical values that these pa- 
rameters may take at low energy. 

The plan of the rest of this paper is as follows. Section 2 introduces the NMSSM 
and the tree-level potential and Higgs boson mass matrices. Section 3 presents an RG 
analysis of the dimensionless couplings between Mqut and MsuSY- In section 4 we 
discuss in general terms the radiative corrections to the Higgs mass matrix using the one- 
loop effective potential. In section 5 we use numerical methods to obtain an upper bound 
on the lightest CP-even neutral Higgs boson mass as a function of the various parameters. 
This is compared to the case in the minimal model. In section 6 we discuss Higgs boson 
masses and couplings in some detail and discuss the effect of radiative corrections. Section 
7 concludes the paper. 

2 The Next-to-Minimal Super symmetric Standard Model 

The most commonly considered supersymmetric model, the MSSM, has, in addition to 
the usual matter and gauge particle content, a Higgs sector containing two Higgs doublet 
superfields Hi and H2. The superpotential is then of the form 

Wmssm = huQu'^H^ + hdQd^Hi + heLe^H^ - fiHiH2, (2.1) 

where generation indices are understood, H1H2 — H^H2 — H2 , with Hj ~ (-f^i, H^), 
H2 = {H2,H2), and the rest of the notation is conventional. In this model the physical 
Higgs spectrum consists of two CP-even and one CP-odd neutral scalars, and the lightest 
neutral scalar h has a mass which is bounded at tree level by < m| [2] . 

In the NMSSM [8,9,10] the particle content of the MSSM is supplemented by a gauge 
singlet superfield, A^. The superpotential is given by 

Wnmssm = KQu^H^ + hdQd^Hi + heLe^Hi + \NH1H2 - ^N^. (2.2) 

The cubic term in N is necessary to avoid a Peccei-Quinn symmetry which would force the 
existence of a light pseudo-Goldstone mode once the symmetry is broken. However there 
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still remains a TLo, symmetry under which all the matter and Higgs superfields $ transform 
as $ ^ q;$ where o? = \ [9] . Note that we have eliminated \x. This can be justified on the 
grounds of naturalness, and its inclusion would only complicate our analysis. The gauge 
singlet field acquires a vacuum expectation value (vev) which plays the role of the mass 
parameter in the MSSM. 

Unlike the case in non-supersymmetric models, radiative corrections do not generate 
large masses of order the cut-off scale (which for SUSY GUTS is essentially the unification 
scale lO^^GcV), although the inclusion of singlets may cause the destabilisation of the 
hierarchy if there are strong couplings to super-heavy particles such as Higgs colour triplets 
[23] . This is however strongly dependent on the structure of the model at the GUT scale, 
and so we shall not discuss it here. Recently it has also been noted that the inclusion of 
non-renormalisable operators suppressed by powers of the Planck mass may lead to the 
introduction of non-logarithmic divergences, which in turn can generate large mass terms 
for the Higgs bosons of the electroweak theory [24] . This effect can only occur in the case 
of a model which has a singlet, and requires that gravity violate the TL-^ symmetry which is 
respected by the renormalisable operators of the theory. The coefficients of such operators 
have not yet been calculated and their size and importance is unclear. 

In our analysis we shall drop all quark and lepton Yukawa couplings apart from that 
of the top quark so that the superpotential reduces to 

W^MSSM ~ htQm' + AiViJiifs - ^iV', (2.3) 

where the superfield = {tL,hL) contains the left-handed top and bottom quarks, 
and contains the charge conjugate of the right-handed top quark. Adopting the usual 
convention of using the same symbols for both component Higgs fields and superfields, the 
fields ffi, H21 and N develop vevs which may be assumed to be of the form 

<ifi>=f'l^V <H2>^{^\, <N>=x = riy, (2.4) 



where ui, 1^2 and x are real, + I'i — ^ — 174 GeV, and tan/? = ^2/1^1- The low 
energy physical spectrum of the Higgs scalars consists of three CP-even neutral states, 
two CP-odd neutral states, and two charged scalars. A third CP-odd state is a Goldstone 



mode which becomes the longitudinal component of the while a further two charged 
degrees of freedom become those of the l^='=s. 

In addition to the potential which can be derived from the superpotential in the usual 
manner, there is a soft supersymmetry breaking potential of the form 

Koft ^htAtQi^H2 XAxNHiH2 - ^A^N^ 

+ mjj^\Hi\^ + m%jH2\'' + m%\N\^ (2-5) 
+ m||rp + m||6^|2 + m^|g|2 

which leads to the full low energy Higgs potential Vq, where 

Vo = ^XiiHi^Hif + \\2{H2^H2f + (As + U){H^^ Hi){H2^ H2) 

- X4\H2^Hif + X5\Nf\Hif + AeliVl'li^al' 

+ \7{N*'^HtH2 + N'^H^HI) + X8\N\^ (2.6) 
+ mi^\Hi\^ + m2^\H2\^ + ms'|iV|^ 

- m4{HiH2N + H.c.) - \m5{N^ + H.c). 

The quartic couplings Aj and the mass parameters must satisfy the following 
boundary conditions at Msusy 

. _ . _ {al + gl) . _ {gl-gl) 
X1-X2- — - — , As - — - — 

A4 = A^ — A5 = Ae = A^, 



2 

A7 = —Xk, Xs = k^, 



(2.7) 



22 22 22 

™1 = > ™2 =^H2 ' "^S = ™iV> 

m4 = XAx, 1715 = kAk 

where gi, §2 are the usual U{1) and SU{2) gauge couplings of the standard model. 

Since we have three minimisation conditions, = and = 0, we may eliminate 
three of the unknown parameters of the theory, which we choose to be mi, m2, and ms, in 
favour of the three vevs. The remaining masses 1714 and are related to the parameters 
Ax and Af. at Msusy as above. Because we know u (from the W mass) the Higgs sector of 
this model is now parametrised in terms of the six parameters A, k, tan/5, r, 777-4, and m^. 
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We shall take the parameters A, k, Ax and to be real, and A, k to be positive, which 
is a sufl&cient condition for the vacuum to conserve CP and leads to a choice of vacuum in 
which all the vevs x, ui, 1/2 are real and positive [9]. 

The range of the parameters is restricted by the condition that the vacuum does not 
break QED in the Higgs sector, which is not automatic in the NMSSM, and is equivalent 
to the condition that > 0, where rric is the mass of the physical charged Higgs H^. 
Another similar problem is that the vacuum may break QCD in the squark sector, but 
this does not occur for sufficiently small At [9, 25]. Slepton vevs will not be discussed here 
since they can be avoided by an appropriate choice of soft parameters. 

The full 10 X 10 mass squared matrix Ai^ for the scalar fields is simple to derive by 
expressing all of the fields in terms of their real scalar and pseudo-scalar parts 

iV = -^(v?3 + «¥'6) (2.8) 
Hi =^{(p7-i(p9) 
Ht = -^{V8 + i<Pio) 



and then using 



When evaluated at the vevs, 



d(fid(fj 



Ml = (2.9) 



(2.10) 



< ^1 >=V2i'i 

< <f2 > = \/2z^2 

< <^3 > = V2x 
<(pi>=0 1,2,3 

Ai^ breaks down to consist entirely of zeros except in one 3x3 block for the CP-even, 
one 3 X 3 for the CP-odd, and two 2x2 blocks for the charged mass matrices (the CP- 
odd matrix and each of the charged matrices have one zero eigenvalue corresponding to a 
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Goldstone mode, and the two non-zero charged eigenvalues are equal). Thus the tree-level 
neutral CP-even (scalar) mass squared symmetric matrix, in the basis {Hi, H2., N}, is 

/ 2Aiz^i 2(A3 + A4)z/iz/2 2X5XU1 

= 2(A3 + A4)i/ii/2 2A2i/| 2X6X1^2 

\ 2A5a;^'i 2Xqxi'2 'iXgx'^ — m^x 

(tan /3[m4a; — Ara;^] —[77140; — Arx^] —■^[777,40; — 2A7a;^] \ 
-[m^x - Xrx"^] cot (3[m4X - Xjx'^] -|-[r?74X - 2A7a;2] (2.11) 
-!f[m^x-2Xrx^] -!^[m^x - 2Xjx^] ^[m^x] ) 

Similarly the tree-level neutral CP-odd (pseudoscalar) mass-squared symmetric matrix, in 
the basis {Hi, H2., N}, is 

(tan/3[r?74a; — Ayx^] [m^x — Xrx"^] ^[m4X + 2X7x'^] \ 

[7714a; - Ayx^j cot (3[m4X - Xrx^] ^[m^x + 2A7a;2] . (2.12) 

^ [m^x + 2X^x'^] ^ [m^x + 2 Ayx^] Smsx + ^ [ttiax - AXjx'^] ) 

Finally the tree- level charged mass-squared matrix, in the basis {Hi, H2}, is 

= (^^^'^ J^^) (7774a; - Xrx' - X4V1V2). (2.13) 

In the limit A, /c — 0, a; — > 00 with Aa; and kx held fixed, the components do not 
mix with the Hi, H2 components in the mass matrices in (2.11)and (2.12). This is just the 
MSSM limit of the NMSSM. 

3 Renormalisation Group Analysis 

Now let us consider which values may be taken by the dimensionless couplings A, k and 
hi- Above some scale at which supersymmetry becomes a good symmetry, and defining 
t = log/U, where /x is the renormalisation scale, the RG equations for these couplings [25] 
are given by 

Sn'^={2X^ + k^ + lhj-lgl-lgl)X 

87r2|^=(3A2 + 3fc2)fc (3.1) 

^dht ^ 1 2 3^2 _ 8 2 _ 3 2 _ 13 

dt ^2 * 3^^ 2^^ 18^^^ 

where gs is the QCD coupling. Following the analysis of reference [26] , they may be written 

in the suggestive form 



(3.2) 



These have three fixed points in the gaugeless {gi = 0, Vz) hmit: 




Of these fixed points, only the last is infra-red stable. This is illustrated in Figure la, where 
the paths traced out by a number of equally spaced points in the X/ht — k/ht plane with 
gi = and /i^ = 10 at the GUT scale are plotted as the energy scale runs from the GUT 
scale of lO^^GeV to lO^GeV. A point in this plane will flow rapidly towards the central 
valley which passes through all three fixed points, and then more slowly along it to the 
stable fixed point. From this figure, it is clear that, regardless of the high energy values, 
the low energy values of the couplings are likely to be somewhere in this one dimensional 
region, but they may well be nearer an unstable than a stable fixed point. 

So far, however, we have neglected the effects of the large QCD coupling in equation 
(3.2); this may be important as discussed in [26]. When we include the effect of gauge 
couplings in our numerical analysis, a different picture emerges. Although at high energy 
the non-zero QCD coupling is small enough not to make much difference, and the flow is 
much as before, at lower energies this contribution dominates, and the effect is that points 
flow towards and along the valley for a short distance before being forced towards the 
origin at lower energies. However, it is noticeable that the region of parameter space to 
which the couplings flow in this figure is still one dimensional. This is illustrated in Figure 
lb, which again has ht = 10 at the GUT scale. 

Both of these plots have been done in the case where the couplings A and k are large; if 
they are small (< 1 at high energy), then the flow towards the flxed point in the gaugeless 
limit is much slower, and with gauge couplings included most points flow rapidly to the 
origin remaining within a region bounded by the one dimensional valley described above 
and the x- and y-axes, as shown in Figure Ic, where /it = 1 at the GUT scale. Finally, if 
A, k are greater than ht at high energies, then they flow very rapidly towards the origin 
and end up on the boundary of this region. 
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In summary, for fixed ht the renormalisation group flows give a curve in parameter 
space which can be attained consistent with triviahty from the assumption of large cou- 
pUngs at the GUT scale, but any point between this curve and the origin is also possible 
given appropriate initial values of the couplings. 

This analysis, with the assumption that the high energy behaviour of the theory is 
described by a GUT model, or at least that perturbative physics continues up to the 
unification scale of order 10-^^GeV, means that we can find maximum values of A and k 
consistent with a given value of ht- If we take k = (in order to maximise the low energy 
value of A), it is possible to solve for Amax (the value of A at which it just becomes non- 
perturbative at the unification scale) as a function of ht{lTeV) (or equivalently of rrit and 
tan/?). This has been done [11,12], and the results are reproduced in Figure 2. 

4 Radiative Corrections 

There have been a number of attempts to calculate radiative corrections to the tree-level 
results presented in section 2 [13, 14, 15, 16, 17, 18, 19, 20]. In the MSSM, the dominant 
radiative corrections are those resulting from loops of top quarks and stop squarks. These 
corrections have been calculated using a range of diagrammatic, RG and one-loop effective 
potential techniques. Radiative corrections due to the top quark and stop squark for the 
upper 2x2 part of the mass matrix are identical to those of the MSSM if we replace 
fj, by Aa;, and these have been studied in a simple approximation [14]. The full analytic 
corrections to other components of the scalar and pseudo-scalar Higgs mass matrices have 
also been calculated [15, 19, 20]. 

Logarithmic effects from loops of light Higgs bosons have been studied using a full 
RG analysis below a hard SUSY breaking scale [16, 18]. However these analyses may be 
criticised on two counts. Firstly, the assumed spectrum of degenerate sparticles at some 
SUSY scale and Higgs bosons degenerate with the top quark at some lower scale is over- 
simplified. Secondly, the RG analysis fails to pick up non-logarithmic radiative corrections 
depending on soft trilinear couplings. These corrections are known to be significant in 
the squark sector, and the corresponding corrections from the Higgs sector may also be 
important. It is possible to complicate the RG approach to take account of some of these 

10 



effects, by systematically decoupling heavy particles below their mass thresholds, and 
introducing finite shifts in boundary conditions, but the elegance of the RG approach is 
then lost [7]. For the cases where the effects of finite diagrams are of interest, the simplest 
approach is to perform a one-loop effective potential calculation. In reference [19] we 
considered a hybrid approach in which squark corrections to all Higgs boson masses were 
calculated in the framework of the one-loop effective potential, and these corrections were 
grafted on to our previous RG analysis [18] involving light Higgs bosons and a degenerate 
top quark. 

In this section we shall abandon the RG approach completely and perform a full 
one-loop effective potential analysis of radiative corrections involving the top quark, stop 
squarks, Higgs bosons, and Higgsinos. We shall introduce the use of the effective potential 
for calculating radiative corrections to mass matrices, present corrections due to the top 
quark and stop squarks, and discuss the calculation of Higgs and Higgsino corrections. 

4-1 One-Loop Effective Potential 

The full one-loop corrected scalar potential is given by 

Vi=Vo + AVu (4.1) 

where Vq is the tree level potential, and AVi is 

Al'. = ^St.A.^(log^--). (4.2) 

jj, is the MS renormalisation scale, and the supertrace is a trace over all fields which couple 
through the mass matrix and includes a factor (— 1)^-^(27 -|- 1) so that a Weyl fermion 
acquires a factor -2, a real scalar a factor 1, and we must remember appropriate colour 
and flavour factors. Ai^ is the field dependent mass-squared matrix, in which the fields 
left after differentiating Vq to obtain J\4^ are not replaced by their vevs. 

The one-loop minimisation conditions are, of course, = 0, z = 1,2,3. After 
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replacing the fields by their vevs, we have explicitly 

2m\i'i + 2Ai^'^ + 2(A3 + \i)i'iV2 + 2\^x^i'i + 2\'jx^i'2 — m4Xi'2 + V2 ^ 
2m\v2 + 2X2V2 + 2(A3 + \/i)viV2 + 2XqV2X^ + 2X'jVix^ — m^vix + \/2 
2m\x + 2Xr,v\x + 2XqV2^ + AX'jViV2X + 4A8a;^ — 2m4ViV2 — 2m^x^ + \/2 



d(fi 
dAVi 



d(p2 
dAVi 



d(p3 



=0 

vevs 

=0 

vevs 

=0 

vevs 



(4.3) 

from which we can see that 



2 _ 21 ^ 9AVi 



(4.4) 



where m||TL are the values of in terms of the vevs as evaluated at tree-level. Having 
thus supplied the vevs, we extremise the potential by choosing the soft masses to satisfy 
equation (4.4). This does not guarantee that this extremum is a minimum of the potential; 
this can only be done by calculating the physical Higgs spectrum and ensuring that all the 
masses squared are positive. 

The minimum so constructed need not be the global minimum of the Higgs potential. 
Thus we examine the potential containing these (radiatively corrected) mass parameters 
m| to determine whether alternative minima with zero vevs for one or more of the scalar 
fields (fi, ip2, and ips exist. In fact the case where only one of these is zero does not occur 
for any of our choices of parameters. We test such possible minima by comparing the value 
of the potential, including top and stop corrections only for simplicity, to discover whether 
these alternative minima are preferred; if so we discard this point in parameter space. 

The radiatively corrected matrices for the CP-even, CP-odd and charged Higgs scalars 
may be calculated in a straightforward manner. In the basis defined by (2.8), the correction 
to the mass-squared matrix is 

d(/),d(/)j 

In fact, this is only approximately true due to Higgs self-energy corrections; these are 
expected to be small for the lightest Higgs bosons [5]. However we should recall that 
the formula for given in equation (2.9) contains implicit dependence on the mass 
parameters , and so will involve dependence on the first derivative of AV. Thus for 
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example the radiative corrections to the scalar mass-squared matrix whose tree-level 
value is given in equation (2.11) are of form 



(4.6) 



The last term represents the shift in the mass matrix due to the radiative corrections to 
the minimisation conditions. 

Jt..2 Squark and Top Corrections 

The corrections to the mass matrices due to top quark and stop squark loops have been 
calculated elsewhere [14, 15, 19, 20], so we simply give the results, in the notation of [19], 
without further explanation. The field dependent top quark and stop and sbottom squark 
mass matrices are given by 

Ml^ = hU\H',\' + \H+\^) (4.7) 
/ m^ + /i2|if0p htXNH^ + htAtH^ -hjH^H+ \ 



htXNH^ + hjAtH^ m| + /i|(|if^|2 + |iy+|2) htXNH^ - htAtH+ 
-hlH^H^ htXNH^ -htAtH} m?Q + hl\H^\'^ 

V m|/ 



B 

(4.8) 

in the basis {tL^fj^, bz,, b'^}, and a bar denotes complex conjugation. The sbottom squarks 
appear because they contribute to the charged Higgs mass. One eigenvalue of (4.8) is field 
independent and may be discarded. Notice that the form of this equation implies that 

I 9 9 I 

At + Xx cot (3 < ^ (4.9) 

Here m^^ and m^^ are the stop squark mass eigenvalues. 

After considerable algebra we obtain the following corrections to the Higgs mass ma- 
trices. The correction to the CP-odd mass-squared matrix in equation (2.12) is given 
by 5M^, where 



■ tan /3 1 



sin/3 



with 



5M'' = I 1 cot (3 ^ I A^ (4.10) 

sin/3 cos /3 sin /3 cos /3 



^h1{Xx)Atf{ml,mll (4.11) 
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and the function / is defined by 

2 2 



The correction to the CP-even mass-squared matrix in equation (2.11) is given by 
SM^, where 



A?. A?.\ /tan/? -1 



sin f3 



5M^= I A?2 Ai' AI3 + -1 cot/9 -4^ I A^ (4.13) 

Ai3 A23 A33y y sin^ _cos^ sin ^ cos /3 



where the Ajj are given by 



a2 3 4 2/. .aMt + ^^cot/3\2 2 2^ 

^2 ^1 



12 tl tl 

At(Ai + Xx cot /3) 



2 



+ 



2 2 
m| — m| 

12 tl 



' /' 2 2 \\ 



a2 3 4 w, ^f At + Ax COt/j ^2 

t2 tl 



3 

2 ™2 „ . / , , , . ^^ ™2 



h'ivl{\x){\v,)f{r4.ml) (4-14) 



Aoo = — ^/t: 



3 ^4 2A 2At(Ai + Axcot/3) . 
^^zy|(log( — + 2 2 ^log( 



m 



^22 — ■^"^"■f 2 v^b\ 4 — ; T 2 2 "^^K — 2~ 

Stt^ v ^ m? ^ mi — mi ^mi 

12 II II 

3 ,4 2Mt(^t + A2;cot/3)\2 2 2n 



12 tl 



2 _COS^ 2 
^23 — ~ ^12 



^33 



ht^2{^^l) [ ^2 _^2 ) ^Ki'^t,) 



Stt^ * ^ \ mi -mi 

12 tl 



and the function g is defined by 

2 

™y ^ ^2 1^^2 ((< + O l°g(^) + 2< - (4-15) 

t2 tl tl 

Finally, the correction to the charged mass-squared matrix in equation (2.13) is SM^, 
where 
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and 



_(_ \ vevs 



/ifi/|(Aa;)2cot/9 



(m? — m| ) (mf — m? ) (m? — m? ) 

51 12 II 01 tl 12 



/i^z^K Ax)2 cot /3 



(m? — m? ) (mf — m? ) ~'~ (m? — mf ) 

12 Oi 12 tl ' II 12 

/i|z/|(Aa:)2cot/3 
(m? — m? )(m? — mf ) 

^ Ol II ' ^ Ol 12 ' 



/i^(Ax)At 



Ol 



dH^dH^ 



_l_ I vevs 



(4.17) 



(4.18) 



.^.^ Higgs and Higgsino Corrections 

Calculation of the corrections due to Higgs loops is much more cumbersome because the 
field dependent Higgs mass-squared matrix in equation (2.9), used to calculate AVi in 
equation (4.2), is 10 x 10. Since it is impossible to determine the eigenvalues of general 
matrices larger than 4x4 analytically, resort to numerical techniques is essential. As is 
clear from (4.6), we are interested in the first and second derivatives of AVi evaluated at 
the vevs, and these are given by 
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where {m^} is the set of the 10 eigenvalues of Mfj, and we implicitly sum over a. It 
is straightforward numerically to obtain m1\yevs from Aifjlvevs, and since the first and 
second derivatives of evaluated at the vevs may also be obtained numerically, the 
problem reduces to a routine numerical task. 

Similarly, since the Higgsino mass matrix is 5 x 5, Higgsino loop contributions must be 
calculated numerically as above. We have neglected the mixing with gauginos for simplicity, 
as in this way we can eliminate two soft mass parameters. Moreover, if we include gaugino 
corrections, then we should include all gauge boson corrections for consistency, particularly 
since there are typically cancellations between particles and their superpartners, as we shall 
see later. The field dependent Higgsino mass term, with Higginos written as 4-component 
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Majorana spinors, in the basis ip = {H^, H^, N, , , is given by ^ipMfjip, with 
defined by 

Mfj = ReM' - i75lmM', (4.21) 
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The contribution to AVi from Higgsinos is then given by 
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(4.23) 



where the trace is over Dirac as well as internal indices. Note that there is a relative factor 
of 2 between this result and that of [22] because we are working with Majorana spinors. 

The numerical calculations are of course very computationally intensive; however, 
numerical inaccuracies should generally be less than about IGeV, particularly in regions 
of parameter space where all of the particles are light. It appears unlikely that the inclusion 
of two-loop efi^ects or removing any of the approximations in the calculations would give 
corrections significant relative to our uncertainty in the parameters. 

5 The Bound on the Lightest CP-even Higgs 

An upper bound on the lightest neutral CP-even scalar in the NMSSM may be obtained 
from the real symmetric 3x3 neutral scalar mass squared matrix by using the fact that 
its minimum eigenvalue is less than or equal to the minimum eigenvalue of its upper 2x2 
block with the result [10] 



ruh^ < Mz^ + {X^v'^ - Mz') sin^ 2/3 



(5.1) 



The upper bound on nth is determined by the maximum value of X{Msusy), henceforth 
denoted X^ax- As discussed in section 3, the value of A^^a; is obtained by solving the 
SUSY RG equations for the Yukawa couplings /it, A and k in the region Msusy = 1 
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TeV to Mgut = 10^^ GeV [11,26]. There we found that for htiMgusr) = 0.5 - 1.0, 
Xmax = 0.87 - 0.70 and for htiMsusr) ^ 1-06,* A 

max — *■ (with k — always). 

We now wish to evaluate the bound including all the radiative corrections. The 
parameters on which the radiatively corrected bound depends are mt, sin/3, r. A, k, 777.4, 
7775, At, rrif^, and m^^ (sbottom masses do not affect the bound, because they only occur 
in the radiative corrections to the charged Higgs mass as discussed in reference [19]). We 
adopt the following approach to these variables. A will, for given rrit and sin/3, be given 
by its maximum value consistent with remaining perturbative, A^ax- Strictly speaking k 
and 7775 should then be zero, but since this results in an axion and thus an unacceptable 
spectrum, we shall instead use the value 0.1 for k. In fact sensitivity of the bound (and 
of Amax) to k and is negligible when they are both small (0.1 is sufficiently small for 
k, while small for means of order tens of GeV). This should however be borne in mind 
when we come to discuss the spectrum, where the lightest scalar is quite sensitive to these 
variables. Of the remaining parameters, 777t is taken as an input parameter, as is the 
heavier stop mass rrif^ . 

The bound also depends on r which may be saturated by requiring that the larger 
eigenvalue of the upper 2x2 block is also an eigenvalue of the whole 3x3 matrix. This 
can be done by performing a unitary transformation consisting of an upper 2x2 unitary 
matrix, parametrised by an angle 9, together with 33 element equal to unity and all other 
elements equal to zero. The condition that the 13 element of the transformed matrix be 
zero leads to the constraint 

cos ^Mfg + sin ^M|3 = (5.2) 

where Mf^ are components of the mass matrix in (2.11). While this guarantees that 
the bound must be saturated, it is possible that there is another lighter mass eigenstate 
consisting of a mixture of singlet and doublet states. This state will certainly depend 
strongly on the values of k and 7775. Examples of this will be discussed in section 6. 

We have thus succeeded in reducing our parameter space to four variables, sin/3, rrif^, 

* ht^Msusv) < 1-06 is the triviality bound, which, together with mt = ht{mt)v sin P, 
where ht{mt) < 1.12, implies the bound mt < 195 GeV. 
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At, and 7714. These variables will all be determined numerically using the Nelder-Mead 
simplex method [27] to maximise the bound as a function of the parameters, imposing the 
constraint (4.9). 

At this point it is worth recalling that in the MSSM the relevant parameters may be 
taken to be tan/3 and rUc- Using the fact that in the minimal model the bound is always 
maximised for sin/3 = 1, the bound is then equal to the 22 element of the mass matrix 

+ SM'^ in (2.11) and (4.13), and is given by 



where the function g is defined in (4.15). Since the coefficient for the Af term is always 
negative * the bound will be maximised for some non-zero value of A^. The constraint 
(4.9) must also be implemented so that value of At will not always reach the value which 
maximises the polynomial. Finally, we note that if we allow m^^ to become comparable 
in size to the top mass, then the coefiicient of the second logarithm in equation (5.3) will 
become negative, and so the bound will be maximised when the two squarks are degenerate; 
this typically occurs only for very large top mass, even if we allow mr to be as low as 



Let us now proceed to the discussion of our results for the bound in the NMSSM (and 
MSSM for comparison) with all radiative corrections included. Figure 3 shows the values 
of sin/3, ht, and A for which the bound is saturated as a function of the top mass, derived 
as described above. Although this figure is for m^^ = 500GeV, the results which can be 
derived for other values of m^^ are virtually identical. This figure is very similar to Figure 
1 of reference [18], which was, however, derived without the inclusion of squark effects, as 
discussed there. 

* Note that ^ < and g/{m'~^ — rn'-^'^ remains finite as m^^ tends to mi^. 




(5.3) 




)g{mt^,mt^) 



250GeV. 
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Values of the bound as a function of the top mass with the heavier stop mass m^^ 
taken to be 500 and 1000 GeV are shown in Figure 4. It is clear from this graph that 
for low top mass the dependence on the stop mass is very slight (as should be obvious, 
as here ht is small), while for larger top masses this dependence is far greater, as the 
radiative corrections are very large. At tree level the decreasing value of A means that 
the bound is a monotonically decreasing function of mt, but for large enough m^^ ^^e 
bound actually starts increasing again with increasing top mass (both the top corrections 



is large). Numerically the bound is around 150 to 155GeV, and is maximised for small 
top mass; but for top masses near the triviality bound and stop masses up to ITeV this 
same value is also approached. The radiative corrections to the bound from Higgs and 
Higgsino loops are relatively small, of order a few GeV, and since they are of opposite sign 
(including Higgs loops decreases the bound, while including Higgsinos increases it) the 
total effect is typically only one or two GeV. In general the optimum value of At is of order 
1.6TeV to 2TeV for m^-^ = ITeV and of order 600GeV to SOOGeV for mf^ = 500GeV. The 
dependence of this optimum value on the top mass is very small, but generally it decreases 
as rrit increases. Typically the value of m^^ which maximises the bound is around a third 
to three quarters the value of m^^ . 

The MSSM bound for the same values of the stop masses is also shown in Figure 4 
for comparison. In this figure Higgs and Higgsino radiative corrections are included in 
addition to the top quark and stop squark corrections in (5.3). Note that, unlike the case 
in the NMSSM, the MSSM bound is a monotonically increasing function of top mass, and 
that it approaches the NMSSM bound for sufficiently large m^. 

The bound is also presented in Figures 5 but here as a function of rrif^, the lighter 
stop mass, for mt = 150GeV and m^^ =500GeV and lOOOGeV respectively. As can be 
seen from these graphs, there is a gradual rise in the bound until it reaches the maximum 
before falling off more rapidly. The very rapid drop in the bound as approaches m^^ is 
caused by the constraint (4.9), which forces + Xx cot (3 to zero in this region. The most 
interesting feature of these figures is that they suggest that an arbitrarily chosen value of 
the lighter stop mass is not unlikely to give a value of the bound close to the maximum 




dependence and the stop corrections can become large only if m^. 
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possible, since the variation of the bound with m^^ is relatively small over much of its 
range. 

6 Higgs Boson Phenomenology 

6.1 Preliminary Remarks 

In this section we shall explore some of the phenomenological implications of our results. 
The qualitative features of certain aspects of our discussion are well known at the tree-level 
[9] ; however we are now in a position to re-examine some of these issues in the light of our 
treatment of radiative corrections. 

We have already emphasised that the lightest CP-even neutral scalar h (the analogue 
of the Higgs boson of the standard model) may be heavier in the NMSSM than in the 
MSSM, although as mt — > 190GeV the two bounds coalesce as shown in Figure 4. However 
from a pragmatic standpoint the key question is whether m/j can exceed Mz and so 
evade discovery at LEPII. This can happen in both the MSSM and the NMSSM due 
to radiative corrections and thus failure to observe h at LEPII cannot be interpreted as 
evidence for favouring the NMSSM over the MSSM. Also we observe that the NMSSM 
bound nih < 150GeV means that if h is not found at LEPII it must necessarily be in 
the intermediate mass region Mz — 2Mw {i.e., too heavy to be produced at LEPII but 
too light to decay into pairs of Ws or Zs) in both models. Thus the upper bound on rrih 
may not be particularly helpful in enabling us to distinguish between the MSSM and the 
NMSSM. 

We can turn the argument of the preceding paragraph around and ask whether the 
observation of h would enable the MSSM and NMSSM to be distinguished. This question 
is analogous to that of distinguishing between the minimal standard model (MSM) and 
the MSSM from an observation of the Higgs boson, and has been widely considered [2,28]. 
The answer depends in part on measuring the ZZh coupling which contains a factor of 
Rzzh — sin{(3 — a) relative to the usual standard model coupling, where a is a mixing 
angle which results from diagonalising the scalar 2x2 matrix. Although Rzzh may be 
small, the ZAh coupling contains a factor of RzAh = cos{P — a) so that both couplings 
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cannot simultaneously be small. In the NMSSM the situation is not so simple since the 
ZZh coupling is derived from diagonalising the scalar 3x3 matrix and so Rzzh is more 
complicated [9]. Similarly the ZAh coupling, where A is the lightest pseudoscalar, in the 
NMSSM will contain a more complicated factor Rzah [9]. 

Given that in the NMSSM physical Higgs boson spectrum there are three CP-even 
scalars and two CP-odd pseudoscalars, one more in each case than in the MSSM, it may 
seem at first sight that the Higgses of the NMSSM, being more numerous, are therefore 
easier to discover. Unfortunately this is not so. The only reason why there are more 
neutral states is due to the complex singlet N , whose scalar and pseudoscalar components 
mix with the neutral components of H\ and H2. Since N has no gauge couplings, this has 
the effect of diluting the couplings of the neutral Higgs particles: there are more of them 
but they all couple more weakly, making them harder to produce. However the dilution 
of the neutral Higgs couplings may be the key to distinguishing between the NMSSM and 
the MSSM, since in the NMSSM there is the possibility of a very light Higgs boson waiting 
to be discovered with higher statistics data from LEPI. This is an intriguing possibility, 
unique to the NMSSM, because although it is possible to have Rzzh = sin{f3 — a) in 
the MSSM, this is always accompanied by RzAh = cos{/3 ~ a) Ki 1 and since small rrih is 
associated with small niA the light decoupled Higgs scenario in the MSSM is ruled out. 
Thus the discovery of a light weakly coupled Higgs boson would be exactly the sort of 
signature which would suggest the NMSSM. 

The charged scalars remain unaffected by the presence of the additional singlet, and 
their gauge and matter couplings in the NMSSM are identical to those of the MSSM. How- 
ever in contrast to the MSSM bound rric > (in the absence of radiative corrections), in 
the NMSSM we may have rric < Mw- There will however be a constraint on how light the 
charged scalars can be since they will give a contribution to 6 — > S7 via a penguin diagram 
involving an intermediate top quark line [29] . Such contributions may be partially or com- 
pletely cancelled by other SUSY diagrams [30]. It would therefore be prudent to keep an 
open mind on the existence of light charged scalars in the range nic — 45GeV — My^ which 
would be accessible to LEPII. We emphasise that this range is not consistent with the 
MSSM and so the discovery of charged Higgs bosons at LEPII would be, if not a smoking 
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gun of the NMSSM, then at least one which is loaded and cocked. 

The Higgs boson mass spectrum of the MSSM at tree level may be expressed in terms 
of the parameter set {mA,tan/3}, so that at tree-level its properties may be expressed as 
contour plots in the uia, tan/3 plane [28]. In the NMSSM, the Higgs boson mass spectrum 
may be expressed in terms of the parameter set {X,k,mc, Ak,tanP,r}. These masses are 
conventionally plotted as a function of rric for a particular choice of A, k, Aj~, tan r. There 
are thus four additional parameters in the NMSSM which may be taken to be A, k, and r 
which are simply not present in the MSSM. For these four additional parameters, reference 
[9] uses A = 0.87, k = 0.63 corresponding to an approximate fixed point, and varies A^ 
over allowed ranges to produce bands of Higgs boson masses, for three difi^erent r values 
of 0.1,1,10. Although it is clear from Figure la that there is indeed a fixed point in the 
gaugeless limit, it is equally clear from Figures lb and Ic that with the QCD and other 
gauge couplings switched on this fixed point is washed away. Instead there seems to be a 
fringe of preferred values of A, k roughly corresponding to a quarter-circle in the X-k plane. 

If one is tempted to despair at the six dimensional parameter space of the NMSSM 
compared to the two dimensional parameter space of the MSSM, it is worth remembering 
that when radiative corrections are taken into consideration things get much worse. In 
addition to ruA and tan/3, the MSSM Higgs sector depends on five additional parameters 
rrit, fi, rrif^, rrif^ and A-t, even neglecting the bottom quark Yukawa coupling. Actually the 
NMSSM fares slightly better since with radiative corrections due to top quarks and squarks 
included there are only four additional parameters mt, m^^, m^^ and At since there is no 
parameter. The Higgs and Higgsino radiative corrections do not introduce any further 
parameters since their tree-level masses, which are used to generate the one-loop effective 
potential, are determined in terms of the tree-level parameters. This would of course not 
be the case if we had included the full effects of Higgsino-gaugino mixing, introducing 
extra gaugino mass paramters. One way to reduce the number of parameters is to assume 
universal soft parameters at the GUT scale and this has been done in the NMSSM [9, 31]. 
Here we shall not restrict the soft parameters in this way, and instead attempt a crude 
exploration of the full parameter space without any constraints from the GUT scale apart 
from the usual SUSY desert assumption that no couplings blow up below Mqut- 
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We shall plot Higgs boson masses as a function of rric over the regions of rric which 
correspond to all Higgs boson masses being greater than zero at tree-level. This turns out 
to be a non-trivial restriction, usually (but not always) due to the lightest scalar mass nih 
diving to zero at the edges of the region. Although radiative corrections tend to expand 
this region, since typically they serve to increase Higgs boson masses, we use the tree-level 
restriction since it is the tree-level mass matrices which are used in the calculation of the 
one-loop effective potential. Charge breaking in the Higgs sector is trivial to check for 
since it just corresponds to becoming negative. The possibility of squark and slepton 
vevs as discussed earlier will not concern us here. However we search for alternative vacua 
in which one or more of vi,V2,x take zero values corresponding to a deeper minimum of 
the (radiatively corrected) potential, as discussed in section 4. This can be an important 
constraint on the allowed range of rric- 

6.2 Higgs Boson Masses and Couplings 

With the preliminaries over we now embark on our exploration of the ten dimensional 
parameter space which characterises the Higgs boson spectrum of the NMSSM. Clearly a 
full survey of Higgs boson masses over all of parameter space is not feasible. The best one 
can do is to define a set of baseline parameters and explore the effect of varying each of 
the parameters in turn. Our baseline parameters are 

A = 0.65 k = 0.1 Ak^O tan/3 = 1.7 r = l 

(6.1) 

mt = ISOGeV m^-^ = 150GeV m^-^ = SOOGeV At = 

We choose these particular values of A, k, tan/3 because, for a top mass of 150GeV, they 
are approximately those which saturate the bound on nih as discussed in section 5, and are 
therefore associated with heavier values of the lightest CP-even scalar. Of the remaining 
parameters, Ak = At = 0, r = 1, are chosen for simplicity, and mt = 150GeV is selected 
as a typical value. The stop masses (with m^^ chosen by convention to be the heavier) are 
again familiar from our discussion of the bound, and ensure that the condition on the stop 
eigenvalues in equation (4.9) is always satisfied over our ranges of parameters, with the 
top quark mass equal in mass to the lightest stop. In all our plots we shall take Ai^ = 0, in 
accordance with the baseline parameter set defined above. In fact, there is often very little 
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freedom of choice of Af. since it is clear from the tree-level mass matrices (2.11) and (2.12), 
that allowing a large positive value of (recall that = kAk) forces the scalar mass 
matrix to have a negative eigenvalue, while allowing a large negative eigenvalue forces the 
pseudo-scalar mass matrix to have a negative value. This is further compounded by the 
structure of the potential, which means that for large positive values of the coefficient 
of the term is large and negative, and so the alternate minimum with x ^ and 
1^1 = 1^2 = will be preferred. 

Figure 6a shows the Higgs boson masses for the set of parameters defined in equation 
(6.1), including radiative corrections due to loops of top quarks, squarks, Higgs bosons 
and Higgsinos. Over the charged Higgs mass range rric = 200 — 250GeV the lightest scalar 
mass (lowest solid) varies between m/^ f=:i 50 — 90GeV while the lightest pseudoscalar mass 
(lower dashed) is approximately 80 GeV. For smaller values of rric than plotted the lightest 
scalar mass dives down to zero, while for larger rric than plotted vacuum with Vi, x = is 
preferred, and these two constraints thus define the allowed region. In Figure 6b we display 
the amplitude of N contained in the lightest scalar (dots), the next-to- lightest scalar (dot- 
dash) and the lightest pseudoscalar (double-dot-dash). It is clear that the lightest scalar 
is predominantly A^, especially near rric = 'iS5GeV, while the next-to-lightest scalar shows 
the opposite behaviour, with a very small component, especially where the lightest 
scalar is predominantly N. The lightest pseudoscalar is nearly entirely singlet. Also shown 
in Figure 6b are the couplings Rzzh (solid) and RzhA (dashed) which confirm that for 
rric ~ 235GeV the lightest scalar h and pseudoscalar A effectively decouple. This behaviour 
was also noted in reference [18] where it was pointed out that the next-to-lightest scalar 
must and does respect the bound at the point where the lightest scalar decouples. However 
since its mass exceeds 100 GeV over the whole region, the next-to-lightest scalar will not 
be visible at LEPII. The lightest scalar may be accessible to LEPI near the left-hand end 
of the region, however. 

As an example of a region of parameter space which is excluded by LEPI, in Figure 7a 
we show the Higgs boson mass spectrum for a value of r = 0.1, with all the other parameters 
set equal to their baseline values as in Figure 6. Smaller r is associated with smaller 
charged Higgs masses, so that the LEPI bound rric > 45GeF is sufficient to exclude the 
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left-hand half of the region immediately. The constraint from b ^ may exclude the 
whole of the region as discussed above, but even without this constraint the right hand 
half of the region is probably excluded from LEPI Higgs searches, as we now discuss. The 
lightest scalar has a mass of around 50 GeV or less, and the lightest pseudoscalar has a 
mass of about 20 GeV, and both these particles have very little singlet component, as 
shown in Figure 7b. The heaviest scalar and the other pseudo-scalar primarily consist 
of singlet, and are approximately degenerate at around half a TeV, playing no important 
phenomeno logical role at LEP. However the LEP complementary searches for Z ^ Z*h 
and Z hA probably eliminate the entire allowed region of masses, since both Rzzh 
(solid) and RzhA (dash) in Figure 7b are quite sizeable on the right-hand part of the 
region. It therefore seems likely that the whole of the region is either excluded or on the 
point of being excluded by LEPI. 

In Figure 8a we show the spectrum of Higgs bosons for r = 5 and the other parameters 
set equal to their baseline values. In this example the lightest scalar h contains very little 
N component particularly in the middle part of the range, according to Figure 8b, and 
couples to the Z almost identically to the Higgs boson of the MSM, according to the solid 
line in the figure. The lightest scalar mass exceeds 100 GeV over much of the region, and so 
is inaccessible to LEPII except at the extreme ends of the range. The lightest pseudoscalar 
and second lightest scalar are both predominantly N . The charged scalar mass, and the 
heaviest scalar and pseudoscalar masses exceed 1 TeV over the entire range for this value 
of r. Thus this large r parameter set produces a Higgs spectrum which will not be easy to 
study even at LHC/SSC. 

All the above plots are for tan/3 = 1.7. In Figure 9a we show the Higgs masses for 
tan/3 = 10 with all the other parameters set equal to their baseline values. Comparing 
Figure 9a to Figure 6a we see that one of the main effects of choosing a larger value of 
tan (3 is substantially to increase the values of rric over which an allowed solution exits, as 
we would expect from the tree-level mass matrix (2.13). Thus we see that larger values 
of rric are associated with larger values of tan/3, as well as with larger values of r as 
previously noted. With rric now in the TeV range the lightest scalar mass is pulled down 
to 30-50 GeV although since it consists almost completely of N and couples weakly to 
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the Z (see Figure 9b), it is not ruled out by current Higgs searches at LEPI. This is 
a good example of a light weakly coupled Higgs boson, characteristic of the NMSSM. 
High statistics LEP data would be required to discover this weakly coupled Higgs boson. 
Similarly the lightest pscudoscalar in Figure 9a, although lighter than its counterpart 
in Figure 6a, remains decoupled and undetectable according to Figure 9b. The physically 
relevant next-to- lightest scalar in Figure 9a remains above 100 GeV, out of range of LEPH. 

Next we turn to the question of the effect of varying k and A. For the above plots 
we selected A; = 0.1 and A = 0.65 because these values served to maximise the lightest 
scalar mass in our analysis of the bound in section 5. Now we shall explore the effect of 
choosing different values of k and A, and in this discussion we shall be guided by our study 
of RG flows in section 3. By comparing Figures la and lb it is clear that the fixed point 
A = 0.87ht and k = 0.63ht is blown away by gauge (in particular QCD) effects. Instead 
we are left with a preferred fringe in the {X/ht) and (k/ht) plane of Figure lb roughly 
corresponding to a circle of radius fa 0.6. This plot is for large ht at the GUT scale, 
corresponding to /j.^ ~ 1.1 at low energy; however Figure Ic shows similar behaviour for 
ht = I at the GUT scale. Our baseline parameters have rrit = 150GeV and tan/? =1.7 
corresponding to ht = 0.92 at a scale of ITeV; this is the value by which the points on the 
fringe of Figure lb must be multiplied to obtain the values of A and k. 

Following the discussion of the preceeding paragraph, a central point on the fringe of 
(X/ht) = {k/ht) = 0-44 corresponds to A = A; = 0.4 and the spectrum with these values 
for the couplings and all other parameters set equal to their baseline values is plotted in 
Figure 10a. Compared to the baseline values of A = 0.65, k = 0.1 in Figure 6a, Figure 10a 
shows nic values of 120-200 GeV, smaller by about lOOGeV, with the left hand end of the 
range determined by the lightest pseudoscalar mass dropping to zero, and the right hand 
end by the lightest scalar mass dropping to zero. The latter has a mass of up to 80 GeV, 
and Figure 10b shows that the lightest scalar (dots) contains little singlet resulting in a 
SM-type ZZh coupling (solid) in contrast to Figure 6b. This puts the lightest scalar within 
the LEP range, with an excluded region corresponding to the LEP limit m/i < 60GeV 
advancing to the left. As before, the lightest pseudoscalar is mainly singlet and so has 
small physical couplings (dashes in Figure 10b). The heavier Higgses are somewhat lighter 
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than before but out of range of LEP. 

Figure 11a shows the Higgs mass spectrum for A = 0.1 and k = 0.6, at the other 
end of the fringe to the basehnc values of Figure 6a. Note that Figure 11a is similar to 
Figure 10a, but the lightest scalar is lighter at around 30 GeV, and since it contains very- 
little N and couples in a standard Higgs like way to the Z (Figure lib) it is excluded 
by present Higgs searches. We conclude that the small A, large k region is not preferred 
phenomenologically. 

Having explored all the regions of the fringe of preferred values of A and k in Figure lb, 
we shall now consider taking these parameters both to be very small. The RG equations 
in section 3 show that A = A; = is an acceptable choice, since if A and k are small at 
the GUT scale, they may still be small at low energies. Furthermore, it is known that in 
the limit A ^ 0, — > 0, r ^ cx), with Ar and kr held fixed, the NMSSM Higgs sector 
reduces to that of the MSSM, so this region is of intrinsic interest for comparisons to the 
minimal model. In Figure 12a we show the Higgs spectrum of the NMSSM for A = 0.1, 
A; = 0.1, r = 5 with all other parameters set equal to their baseline values. This spectrum 
should and does resemble that of the MSSM. For one thing the charged mass exceeds 
MiY, as it should. Also it is clear from Figure 12b that the lightest pseudoscalar (double- 
dot-dash) and the second lightest scalar (single-dot-dash) can be identified with to a 
good approximation over most of the region. The remaining states of Figure 12a then 
correspond to the MSSM states. The lightest scalar of mass 50 GeV couples in the SM 
way (solid curve in Figure 12b) and so is already excluded by LEPI. Additional radiative 
corrections due to non-zero At values and/or larger top and stop masses could increase 
the lightest scalar mass beyond current limits. Such additional radiative corrections are 
discussed in section 6.3. 

We can take the MSSM limit of our results in a much more direct way by removing 
the N components of the Higgs mass matrices by hand, and numerically taking A and k 
to be very small, and r to be very large, with fi selected arbitrarily and appearing in the 
radiative corrections to the 2x2 mass matrices instead of Xx. The resulting spectrum 
in Figure 13a, for = and tan/3 = 1.7 is quite similar to that in Figure 12a, once the 
lightest pseudoscalar and second lightest scalar of that plot have been removed. The range 
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of rric is extended with rric > 80 GeV as expected. The hghtest scalar again has a mass 
ruh ~ 50GeV with Rzzh = sin(/3 — a) 1. 

In Figure 13b we show the corresponding MSSM plots with tan (3 = 10, and all other 
parameters as before. The limit of tan/3 ^ oo is of course where the MSSM bound is 
saturated. In this limit the 2x2 scalar mass-squared matrix is diagonal with elements m\ 
and M| (plus radiative corrections). Furthermore, = m\ + at tree-level. This 
explains why the spectrum in Figure 13b consists of a flat scalar mass plus a rising scalar 
mass degenerate with the pseudoscalar. In this case, since /3 ~ 7r/2 and o: ~ (as the scalar 
matrix is diagonal), it is clear that Rzzh = sin{j3 — ck) ~ 1 and RzAh = cos{j3 — ck) ~ 0, 
where h is defined to be the flat curve in Figure 13b {id est, for values of rric less than the 
point at which the two scalar curves bounce off each other, the lighter of the two scalars 
is decoupled, while for larger values of mc the converse is true.) 

6.3 The Effect of Radiative Corrections 

All the spectra plotted in section 6.2 include top quark and squark radiative corrections 
as well as Higgs and Higgsino radiative corrections. However in all these plots Ai = 0, 
with rrit = m^^ = IbOGeV, m^^ = bOOGeV. We shall now examine the effect of stop squark 
spectra with At ^ 0. To understand how important the effect of non-zero At may be, we 
shall begin by considering the effects of each of the other radiative corrections in turn. 
Since these effects are far harder to compare accurately on a logarithmic scale, we shall 
restrict our discussion to the lightest Higgs boson masses which we shall plot on a linear 
scale. We shall use the baseline parameters deflned in (6.1), and for comparison we shall 
in each case plot the tree-level masses as dotted lines, so that the effect of the particular 
radiative correction is easily seen. 

Figure 14 shows the lightest Higgs boson masses corresponding to the baseline param- 
eters above, so that the Higgs masses are identical to those in Figure 6a. The corresponding 
tree-level masses without top, stop, Higgs or Higgsino radiative corrections are indicated 
by dotted lines. The lightest scalar mass is unaffected by radiative corrections at its peak, 
where it decouples, but at this point the corrections are largest for the next-to-lightest 
scalar. Away from this point the corrections are shared between the two scalars, with 
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the pseudoscalar - which is only weakly coupled across the whole range - receiving small 
corrections. Radiative corrections in this case increase the second lightest CP-even scalar 
mass by as much as 7 GeV. 

Figure 15 shows the effect of radiative corrections due to top quark loops alone, cor- 
responding to the parameters and notation of Figure 14, but with the Higgs and Higgsino 
radiative corrections subtracted. It is clear by comparing Figure 15 to Figure 14 that the 
Higgs and Higgsino corrections have relatively little effect compared to top quark loops, 
the most noticeable change being felt by the pseudoscalar mass. As we shall see, the Higgs 
and Higgsino radiative corrections are individually quite large but tend to cancel between 
themselves (the cancellation being exact in the limit of exact supersymmetry). Top quark 
loops are therefore responsible for the bulk of the radiative corrections for the baseline 
parameters. 

Figure 16 shows the effect of radiative corrections due to Higgs bosons alone, corre- 
sponding to the parameters and notation of Figure 14, but this time with only the Higgs 
boson radiative corrections included. Compared to the tree-level results (dots), the Higgs 
radiative corrections can be quite substantial and tend to depress the scalar masses and 
enhance the pseudoscalar mass. Figure 17 shows the effect of radiative corrections due to 
Higgsinos alone for the usual parameters and notation. In Figure 17 the comparison to 
the tree-level dotted lines shows that Higgsinos tend to enhance scalar masses in such a 
way as to partially cancel the effects of Higgs corrections in Figure 16. The pseudoscalar 
mass is also enhanced in this case, however. 

We are now ready to examine the radiative corrections due to squark spectra with 
Af 0. Again we use the baseline parameters in (6.1), apart from the following two sets 
of stop parameters which were encountered previously in our discussion of the bound in 
section 5: 

(a) rut^ = 150GeV,mt-^ = 500GeV,^t = 700GeV 

(b) m^^ = 600GeV, m^^ = ITeV, At = l.STeV 

Parameter set (a) uses the same stop masses as in the baseline parameters but now 
involves a non-zero At value chosen to maximise the effect of radiative corrections. Pa- 
rameter set (b) involves heavier stop masses and larger At value, again chosen to maximise 
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the effects of radiative corrections. In Figure 18a we show the Higgs spectrum for param- 
eter set (a), including the radiative corrections due to the squark spectrum as well as top, 
Higgs and Higgsino loops, with the corresponding tree-level spectrum again represented by 
dotted lines. We see that the effect of the At = 700 GeV value in Figure 18a is dramatic 
for the physical second lightest scalar with radiative corrections to its mass of up to 18 
GeV as compared to 7 GeV with = in Figure 14 (with all other parameters the same 
in the two cases). 

In Figure 18b we show the spectrum for the heavier squark spectrum set (b) above, 
again including all other radiative corrections in this plot and comparing it to the baseline 
tree-level masses represented by dots. The physical second lightest scalar in Figure 18b 
has radiative corrections to its mass of up to 25 GeV as compared to 18 GeV in Figure 18a. 
We conclude that the effect of choosing larger stop masses and large non-zero At values 
may significantly raise the mass of the lightest physical Higgs bosons from their values 
shown in the plots of section 6.2. 

7 Conclusion 

We have performed a comprehensive analysis of the radiative corrections to the Higgs 
boson mass spectrum in the NMSSM using the one-loop effective potential. Analytic 
results for top quark and squark corrections are reviewed, and our numerical procedure 
for including Higgs and Higgsino corrections is described. The bound on the lightest CP- 
even scalar, including the radiative corrections mentioned above, was then discussed. In 
order to find the absolute bound as a function of rrit including all radiative corrections, we 
maximised over the parameter space using analytic and numerical techniques. Our final 
results are summarised in Figure 4. 

We have presented a numerical analysis of the RG flows for the dimensionless cou- 
plings, relying on the assumption of perturbative behaviour up to Mqut- We emphasise 
the infra-red fixed points of the gaugeless limit are not the relevant points for the low en- 
ergy behaviour because they are washed away by the large QCD corrections, and instead 
the parameter space is limited to a region of the (X/ht) — (k/ht) plane bounded by the 
axes and an approximate quarter circle, as shown in Figure lb. 
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A general discussion of Higgs boson phenomenology including radiative corrections 
was then given. Since the parameter space of the Higgs sector of the NMSSM is multi- 
dimensional, we selected a baseline set of parameters given in (6.1), and discussed the 
effect of varying each of the parameters in turn. Our selection of k and A was based on 
our study of RG flows described above. We found that larger rric values are associated 
with larger values of tan /? and r. The NMSSM with small k and A and large r was also 
compared to the MSSM (see Figures 12 and 13). Possible characteristic signatures of the 
NMSSM include light charged scalars, and weakly coupled light neutral scalars, and there 
are regions of parameter space where either or both signals is possible (see Figures 7 and 
9). Equally there are other regions of parameter space where neither of these characteristic 
signals is present (see Figures 8 and 10). 

Finally Figures 14-18 show the relative sizes of the various radiative corrections. These 
show that the dominant effect is usually that of the top quark in conjunction with the 
stop sector. Figure 18 emphasises the huge radiative corrections which are possible for 
large At and non-degenerate stop squarks. Higgs and Higgsino corrections are typically 
rather smaller, and tend to cancel due to the (softly broken) supersymmetry, as shown 
in Figures 16 and 17. This cancellation also occurs in the top-stop sector, but is less 
noticeable due to the large soft squark masses which are assumed. 

If Higgs bosons are discovered at LEP or LHC/SSC, then an important question will 
be whether they are associated with the MSM, the MSSM, or some other model. If they 
are not associated with the MSSM, then it is possible that they arise from the NMSSM, 
which has a more general structure. The discussion of radiative corrections to Higgs boson 
masses in the NMSSM and the phenomenological discussion presented in this paper will 
help to decide this question. 
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Figure Captions 



Figure la: Renormalisation group flow of a set of points in the {X/ht)-{k/ht) plane as 
the energy scale is varied from lO-'^^GeV to ITeV. /it(lO^^GeV) = 10 and gauge couplings 
have been set to zero. 

Figure lb: Renormalisation group flow of a set of points in the {X/ht)-{k/ht) plane as 
the energy scale is varied from lO^^GeV to ITeV. ht{lO^^GeV) = 10 and gauge couplings 
have been included. 

Figure Ic: Renormalisation group flow of a set of points in the {X/ht)-{k/ht) plane as 
the energy scale is varied from lO^^GeV to ITeV. ht{10^^GeV) — 1 and gauge couplings 
have been included. 

Figure 2: Amax against /it(lTeV) 

Figure 3: The values of sin/3 (long dashes), A (short dashes), and ht (solid line) for 
which the bound is maximised, for the value m^^ = ITeV. Results for other values of rrif^ 
are very similar. 

Figure 4: The bound on the lightest neutral CP-even scalar against m^. The two solid 
lines are for the NMSSM; the dotted lines for the MSSM. In each case the upper (lower) 
line is for m^^ =lTeV (SOOGeV). 

Figure 5: The bound on the lightest scalar against m^^ expressed in units of m^^. The 
upper (lower) line represents the case where m^^^lTeV (500GeV). 

Figure 6a: Higgs boson masses against nic in the NMSSM. The parameters are r = 1.0, 
tan (3 = 1.7, A = 0.65, k = 0.1, Ak = 0, mt = ISOGeV, m^^ - 150 GeV, ruf^ = 500 GeV, 
and At = 0. Solid lines represent CP-even, and dashed lines CP-odd mass eigenstates. 

Figure 6b: Higgs boson couplings and N-amplitudes corresponding to Fig. 6a. The solid 
line represents the Rzzh coupling, the dashed line the RzhA coupling. The amplitude of 
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singlet field N contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar h, a dot-dashed line for the next-to-lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 

Figure 7a: Higgs boson masses against rric in the NMSSM. The parameters are as in 
Fig. 6a except that r = 0.1. Solid lines represent CP-even, and dashed lines CP-odd mass 
eigenstates. 

Figure 7b: Higgs boson couplings and N-amplitudes corresponding to Fig. 7a. The solid 
line represents the Rzzh coupling, the dashed line the Rzha coupling. The amplitude of 
singlet field N contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar /i, a dot-dashed line for the next-to- lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 

Figure 8a: Higgs boson masses against rric in the NMSSM. The parameters are as in 
Fig.6a except that r = 5.0. Solid lines represent CP-even, and dashed lines CP-odd mass 
eigenstates. 

Figure 8b: Higgs boson couplings and N-amplitudes corresponding to Fig. 8a. The solid 
line represents the Rzzh coupling, the dashed line the Rzha coupling. The amplitude of 
singlet field N contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar /i, a dot-dashed line for the next-to- lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 

Figure 9a: Higgs boson masses against rric in the NMSSM. The parameters are as in 
Fig. 6a except that tan/? = 10. Solid lines represent CP-even, and dashed lines CP-odd 
mass eigenstates. 

Figure 9b: Higgs boson couplings and N-amplitudes corresponding to Fig. 9a. The solid 
line represents the Rzzh coupling, the dashed line the RzhA coupling. The amplitude of 
singlet field contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar /i, a dot-dashed line for the next-to- lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 
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Figure 10a: Higgs boson masses against rric in the NMSSM. The parameters are as 
in Fig.6a except that A = 0.4, k = 0.4. Sohd hnes represent CP-even, and dashed hnes 
CP-odd mass eigenstates. 

Figure 10b: Higgs boson couplings and N-amphtudes corresponding to Fig. 10a. The 
solid line represents the Rzzh coupling, the dashed line the RzhA coupling. The amplitude 
of singlet field N contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar h, a dot-dashed line for the next-to-lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 

Figure 11a: Higgs boson masses against rric in the NMSSM. The parameters are as 
in Fig.Ga except that A = 0.1, A; = 0.6. Solid lines represent CP-even, and dashed lines 
CP-odd mass eigenstates. 

Figure lib: Higgs boson couplings and N-amplitudes corresponding to Fig.lla. The solid 
line represents the Rzzh coupling, the dashed line the RzhA coupling. The amplitude of 
singlet field contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar h, a dot-dashed line for the next-to-lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 

Figure 12a: Higgs boson masses against rric in the NMSSM. The parameters are as 
in Fig.6a except that A = 0.1, k = 0.1. Solid lines represent CP-even, and dashed lines 
CP-odd mass eigenstates. 

Figure 12b: Higgs boson couplings and N-amplitudes corresponding to Fig. 12a. The 
solid line represents the Rzzh coupling, the dashed line the RzhA coupling. The amplitude 
of singlet field N contained in each of the mass eigenstates are indicated by a dotted line for 
the lightest scalar h, a dot-dashed line for the next-to-lightest scalar, and a dot-dot-dashed 
line for the lightest pseudo-scalar A. 

Figure 13a: Higgs boson masses against rric for the MSSM. The parameters are tan (3 = 
1.7, mt = ISOGeV, m^^ = 150 GeV, m^^ = 500 GeV, = and = 0. CP-even and 
CP-odd mass eigenstates are represented by solid and dashed lines respectively. 
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Figure 13b: Higgs boson masses against rric for the MSSM. The parameters are as in 
Fig. 13a except that tan/? = 10. CP-even and CP-odd mass eigenstates are represented by 
sohd and dashed hnes respectively. 

Figure 14: Higgs boson masses against rric in the NMSSM. The parameters are r = 1.0, 
tan/? = 1.7, A = 0.65, k = 0.1, Ak = 0, rrit = ISOGeV, m^^ = 150 GcV, m^^ = 500 GeV, 
and At = 0. Solid lines represent CP-even, and dashed lines CP-odd mass eigenstates 
and include radiative corrections due to loops of top quarks, squarks, Higgs bosons and 
Higgsinos. This plot is as in Fig. 6a except that it is rescaled to show the lighter mass 
eigenstates in greater detail. The dotted lines indicate the tree-level spectrum for the 
same parameters. 

Figure 15: Higgs boson masses against rric in the NMSSM. The parameters are as 
in Fig. 14. Solid lines represent CP-even, and dashed lines CP-odd mass eigenstates and 
include radiative corrections due to loops of top quarks and squarks only. The dotted lines 
indicate the tree-level spectrum for the same parameters. 

Figure 16: Higgs boson masses against rric in the NMSSM. The parameters are as 
in Fig. 14. Solid lines represent CP-even, and dashed lines CP-odd mass eigenstates and 
include radiative corrections due to loops of Higgs bosons only. The dotted lines indicate 
the tree- level spectrum for the same parameters. 

Figure 17: Higgs boson masses against rric in the NMSSM. The parameters are as 
in Fig. 14. Solid lines represent CP-even, and dashed lines CP-odd mass eigenstates and 
include radiative corrections due to loops of Higgsinos only. The dotted lines indicate the 
tree- level spectrum for the same parameters. 

Figure 18a: Higgs boson masses against rric in the NMSSM. The parameters are as in 
Fig. 14 except that At = 700GeV. Solid lines represent CP-even, and dashed lines CP-odd 
mass eigenstates and include radiative corrections due to loops of top quarks, squarks, 
Higgs bosons and Higgsinos. The dotted lines indicate the tree- level spectrum for the 
same parameters. 
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Figure 18b: Higgs boson masses against rric in the NMSSM. The parameters are as in 
Fig. 14 except that m^-^ = 600 GeV, m^-^ = 1.0 TeV, At = 1.8 TeV. Sohd hnes represent 

CP-even, and dashed Hnes CP-odd mass eigenstates and include radiative corrections due 
to loops of top quarks, squarks, Higgs bosons and Higgsinos. The dotted lines indicate the 
tree- level spectrum for the same parameters. 
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